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Abstract 

We propose a method to calculate the expectation values of an operator similar to the 
Wilson loop in the large N limit of field theories. We consider AT = 4 3+1 dimensional 
super- Yang-Mills. The prescription involves calculating the area of a fundamental string 
worldsheet in certain supergravity backgrounds. We also consider the case of coincident 
M-theory fivebranes where one is lead to calculating the area of M-theory two-branes. 
We briefly discuss the computation for 2+1 dimensional super- Yang-Mills with sixteen 
supercharges which is non-conformal. In all these cases we calculate the energy of quark- 
antiquark pair. 
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1. Introduction 



It has been expected for some time that the 't Hooft limit [1|] of large iV gauge theories 
is related to a string theory (see and references therein). In a precise string theory 
was proposed for the 't Hooft limit of M = 4 super- Yang-Mills in 3+1 dimensions, based 
on earlier studies [|J. The 't Hooft limit is defined as the limit of iV — ► oo keeping g Y M^ 
fixed. In this limit we get weakly coupled string theory on AdS$ x S 5 where the radius of 
the five-sphere and the curvature radius of anti-de Sitter are proportional to (gyM-^) 1 ^ in 
string units. There is also a flux of the Ramond-Ramond self dual five-form field strength 
on the five-sphere. The string coupling is g ~ 9ym an d goes to zero in the t' Hooft limit. 
In general we do not know how to solve free string theory on AdS§ x S 5 . However, when 
the gN is large the radius of curvature is large and we can use the string in background 
fields approximation. In it was shown how to calculate conformal dimensions of 

operators and correlators in conformal field theory in terms of supergravity when gN is 
large. In this paper we consider the problem of calculating the expectation values of Wilson 
loop operators. The proposal is that these expectation values correspond to the area of 
a worldsheet whose boundary is the loop in question. We will further consider similar 
observables for the M5-brane theory (the conformal (0,2) six dimensional theory). We also 
discuss Wilson loops in non-conformal theories associated with D-twobranes. 



2. The Wilson loop 

Consider a Yang-Mills theory. The Wilson loop operator is 

W(C) = ^TrPe$c A (2.1) 

where C denotes a closed loop in spacetime and the trace is over the fundamental rep- 
resentation. We will be considering mostly the Euclidean field theory. We can view the 
Wilson loop as the phase factor associated to the propagation of a very massive quark in 
the fundamental representation of the gauge group. A loop that is often considered is a 
rectangle as indicated in figure 1 . From the expectation value of this rectangular Wilson 
loop it is possible to read off the energy of a quark-antiquark pair. Namely, in the limit 
T — ► oo the expectation value of the Wilson loop is 



(W(C)) = A(L)e- TE{L) 
1 



(2.2) 



Figure 1: Contour used to extract the quark-antiquark force from the Wil- 
son loop. The vertical direction indicates Euclidean time and the horizontal 
direction indicates one of the spatial coordinates. This contour lives in four 
dimensional Euclidean space. The parameters 6\p are for later reference. 

where E(L) is the energy of the quark-antiquark pair. 

In order to perform this calculation for the cases of interest it will be necessary to 
introduce massive quarks. To this effect consider breaking U(N + 1) — > U(N) x U(l) 
by giving some expectation value $ to a Higgs field. Then the massive W-bosons have 
a mass proportional to |$| and transform in the fundamental representation of U(N). So 
in the limit |<fr| — > oo they provide the very massive quarks necessary to compute Wilson 
loops in the U(N) theory. Notice that we are interested in physics for energy scales much 
lower than |$| so that that the U(N) theory is effectively decoupled from the U(l) theory. 
Consider the equation of motion for the massive W boson. Extracting the leading time 
dependence as W = e~' l ^\ t W we get an equation for W which to first order in 1/|$| reads 

(do-iAo-ie^^W = (2.3) 

where we have defined 6 1 = jj-. Notice that Aq and Xj are matrices in the adjoint of 
U(N). This implies that if we consider this massive W boson describing a closed loop C 
its interaction with the U(N) gauge field will lead to the insertion of the operator 

W{C) = ^TrPj § ds[A » {a) ^ +eI {s)X ' {a) ^ ] (2.4) 

The difference with ( |2.1|) is the fact that we have an extra coupling to X 1 . The operator 
in (|2.4j ) is determined by the contour C (or a At (s)) as well as a function 9(s) mapping 
each point on the loop to a point on the five-sphere. We are interested in this operator 
because it is the one that naturally arises when we consider the propagation of a massive 
W-boson. The appearance of X might seem surprising at first sight, but it is obvious 
when we remember that a string ending on a p-brane is not only a source of electric field 
but it also carries "scalar" charge for the fields X since it is pulling the brane. In fact this 
coupling is crucial to understand the BPS bound for strings stretching between different 
branes ||. In the calculations below 9(s) will be basically constant. 
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3. Relation to supergravity 



A natural proposal for the expectation value of the Wilson loop is 



(W(C)) ~ e" 



(3.1) 



where, in the large gN approximation, S is the proper area of a fundamental string world- 
sheet which at the boundary of AdS describes the loop C and lies along 9 (s) on S 5 . 
See figure 2. In general we should consider the full partition function of string theory on 
AdS§ x S 5 with the condition that a string worldsheet is ending on the loop C and the 
points 8(s) on S 5 at the boundary of AdS. This is a natural proposal in terms of the iden- 
tification proposed in |7]|| for relating gauge theory observables to calculations on AdS. 
However the right hand side in ( |3.1| ) contains also the contribution from the mass of the 
W-boson and it is therefore infinity Subtracting this contribution we find a finite result 
for the Wilson loop operator 



(W(C)) ~ lim e -( s *~ m) 



(3.2) 



Where £ is the total length of the Wilson loop, measured with the flat Minkowski metric 
appropriate to the gauge theory, and $ is the mass of the W-boson. The equation ( |3.2| ) is 
our final recipe for computing the Wilson loop. This result is not "zig-zag" invariant, in 



the sense of 0, since the operator ( |2.4| ) is not invariant, as opposed to (|2.1| ). 




Figure 2: Proposal to calculate Wilson loop expectation values. We should 
consider the partition function of string theory on AdS§ x S 5 with a string 
worldsheet ending on the contour C on the boundary of AdS. 
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4. Quark anti-quark potential 

In this section we consider the calculation of a rectangular Wilson loop as in figure 
1. We take the angle I (s) = 9$ to be a constant. We consider the limit T — > oo. In this 
limit the problem becomes translational invariant along the T direction. We put the quark 
at x = —L/2 and the anti-quark at x = L/2. Here "quark" means an infinitely massive 
W-boson connecting the N branes with one brane which is far away in the direction 9q. 
The action for the string worldsheet is 



S = ^7 J drda yjdet G MN d a X M dpX N (4.1) 



where G mn is the Euclidean AcLSk x S 5 metric 



ds 2 = a' 



U 2 dU 2 

— -(dt 2 + dx l dx l ) + R 2 -^ + R 2 dQ 2 5 



(4.2) 



where R = [A-KgN) 1 ^ is the radius in string units and U = r/a' has dimensions of energy. 
Notice that the factors of a' cancel out in (|4.1|), as they should. Since we are interested in 
a static configuration we take r = t, a = x so that the action becomes 

S = ^ / dx^(d x U) 2 + UyR* (4.3) 

We need to solve the Euler-Lagrange equations for this action. Since the action does not 
depend on x explicitly the solution satisfies 

U 4 

= constant (4-4) 



+ m/R 



Defining Uq to be the minimum value of U, which by symmetry occurs at x = 0, we find 
that the solution is@ 

R 2 [ u / u ° dy 



x = — (4.5) 



where Uq is determined by the condition 

L R 2 f°° dy R 2 V2n 3 / 2 



2 U J 1 y2^43i U T(1/4Y 



(4.6) 



All integrals below can be calculated in terms of Elliptic of Beta functions B . 



M> U 



Figure 3: (a) Initial configuration corresponding to two W-bosons before we 
turn on their coupling to the U(N) gauge theory, (b) Configuration after we 
consider the coupling to the U(N) gauge theory. This configuration minimizes 
the action. The quark-antiquark energy is given by the difference of the total 
length of the strings in (a) and (b). 

The qualitative form of the solution is shown in figure 3. Notice that the string approaches 
the point x = L/2 quickly for large U 

§ - x ~ ^ , U » U . (4.7) 



Now we compute the total energy of the configuration. If we just plug in the solution 
) in (^73|) , we find that the answer is infinity. However as we said above this infinity is 
simply due to the fact that we are including the mass of the W-boson which corresponds to 
a string stretching all the way to U = oo. We can regularize the expression by integrating 
the energy only up to U max . Subtracting the regularized mass of the W-boson which is 
Umax/ (2tt) we find a finite result! 



E 



' 2tt 



dy 



y 



1 -1 



E = - 



47r 2 (2^ M AQ 



1/2 



(4.f 



3 A convenient way to do the integral is to multiply the whole integral by y x , calculate the 
two terms independently as a function of A and then set A = 0. 



We see that the energy goes as 1/L, a fact which is determined by conformal invariance. 
Notice that the energy goes as (gN) 1 / 2 as opposed to gN which is the perturbative result. 
This indicates some screening of the charges. The above calculation makes sense for all 
distances L when gN is large independently of the value of g, this suggest that one could 
define a magnetic Wilson loop operator which for large gN would be determined in terms of 
classical D-string solutions with prescribed boundary conditions at infinity. In the standard 
't Hooft limit the interaction between Wilson loops is governed by g which goes as 1/N. 

4-1. Case of non- constant angle 
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Figure 4: Configuration of a U(N + 2) gauge theory Higgsed to U(N) x 
U(l)i x {7(1)2. This is a view in the transverse space. We also show two 
massive W-bosons they are characterized by the the angle of the Higgs ex- 
pectation value of the {7(1) factor that they are associated to. 

Now we consider the case where the "angle" of the two quarks is different. This arises 
when we break U(N + 2) — > U(N) x f7(l)i x {7(1)2 by giving expectation values $i,$2 
to the two {7(1) factors as indicated in figure 4. Then the angles are Q{ = $i/|$|. So 
we consider a W-boson described by a string going between the iV branes and the brane 
associated to f7(l)i and a W-boson going between the brane associated to {7(1)2 and the 
N branes as indicated in figure 4. Notice that the orientation of the string determines 
whether we have a quark (transforming in the fundamental of U(N)) or an anti-quark 
(transforming in the antifundamental of U(N)). The potential for this configuration can 
be calculated in terms of the large T limit of the expectation value of the Wilson loop 

— * 

shown in figure 1, with different values of 9 on each vertical line. So we should consider a 
string worldsheet which at x = L/2 goes to U = 00 and to the point 6\ of the five-sphere 
and at x = —L/2 goes to U = 00 and to the point 9 2 of the five-sphere. The action for a 
time independent configuration is 

S=^- [dxyj (d x U) 2 + UHdJ) 2 + m/R* . (4.9) 
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From the symmetries of the problem we see that the string will lie along a great circle of 
the sphere. So if we call 9 the angle along this great cicle we can choose #1,2 = ±A#/2. 
The problem then becomes symmetric around x = 0. We can solve the Euler-Lagrange 



equations as above by using the fact that the Lagrangian (|4.9| ) is independent of x and 
9 so that we have conserved quantities associated to "energy" and "angular momentum" 
(interpreting x as time). Solving these equations we find 



x 



U/U 



dy 



9 =1 



U/U 



i/V(y 2 -l)(2/ 2 + l-P) 
dy 



(4.10) 



and the parameters Uq, I are determined by the conditions 



where 



- =x(u = 00) = ^VT^Phii) 

A9 

— =9(U = 00) = ll 2 (l) , 



(4.11) 



mi 
hii) 



1 



(1 -l 2 )V2^P 
1 



[2-l 2 )E 




n 



1- P 

2- f 



TV 



-.F , 

V2^P \ 2 



1-f 



(4.12) 



/ 2 



with F, E are elliptic integrals of the first and second kind. We can also calculate the 
energy of the system, substracting the mass of the W-bosons and we find 



E 



2U 



2tt 



dy 



y 



v / (y 2_ 1)(2/ 2 + 1 _ /2) 

- (2 ^ MiV)1/2 (i-z 2 ) 3/2 A 2 (0 



1 



(4.13) 



Where / is a function of the angle ( |4.11| ). It is interesting to notice that when A# — > n 
then / — - > 1. Then the solution looks like two straight strings going down to U = and the 
energy ( |4.13| ) goes to zero, as expected since this is a BPS configuration. 



6=9] 



V=<x 



u=o 



u=o 



(a) 



(b) 



Figure 5: Solution when the angles associated to the two W-bosons is differ- 
ent, (a) shows the projection on the x, U plane and (b) shows the projection 
on the U, 9 plane, where U now is the radial distance. 

5. M-theory membranes 

If we study the theory of coincident M-theory fivebranes, the (0,2) conformal field 
theory in six dimensions ||10|| , we are led to consider M-theory on AdSj x S 4 . In this 



case one could define Wilson "surface" observables IfLlf . Since we do not have an explicit 
formulation of the theory, we do not have a formula analogous to ([T4|). However we 
could define the Wilson "surfaces" as the phase factor associated to the propagation of 
a very heavy string on branes (sustracting the part proportional to the free propagation 
of the heavy string). In order to be more precise, let us suppose that we start with 
N + 1 branes and they we Higgs by separating one of the branes. A membrane stretched 
between the iV five-branes at the origin and the Higgsed five-brane behaves as a string with 
tension proportional to the separation of the branes. We could consider this heavy string 
as probe for the unbroken conformal field theory associated with the iV branes that are 
still together. The procedure is analogous to what we saw above. The Wilson "surface" 
operator is defined to be the extra phase factor associated with the interaction of the heavy 
string with the N fivebranes. This Wilson "area" operator in the supergravity picture is 
defined by requiring that a membrane ends at the boundary of AdSj x S 14 on the surface 
that defines the operator. Notice that we also have to specify a map from the surface to 
S 14 for the same reasons described above for M = 4 super- Yang-Mills. Again we substract 
the term corresponding to the free propagation of the heavy string to obtain a finite result. 
For large iV we can trust the supergravity result. 



As an example, consider a pair of parallel, infinite strings corresponding to membranes 
ending on the fivebrane. Let us choose them with opposite orientation but in the same 
direction on S 4 . This problem is translational invariant along time and the direction of 
the strings. So the problem of determining the minimal 3-surface reduces, as above, to 
finding the minimum of the action 

rpj I p 

S = / dxy/(dV)* + V 3 /R 3 (5.1) 

(27IT J 

where now R 3 = nN and V = r/l 3 has dimensions of (energy) 2 . The strings have length 
L' and are separated by a distance L in the direction x. We obtain the solution 

(5.2, 

where 

2 "l/ 1/2 r (D 

If we calculate the energy we find 



(5.4) 



E N 8v^Fr(f) 3 

r(i)3 

The dependence on L is the one expected from conformal invariance. 

5.1. Wilson loops in non- conformal theories 

Consider 2+1 dimensional super- Yang-Mills with sixteen supercharges which is the 
theory describing coincident D2 branes. We can define the Wilson loop operator as in 
( |2.4j) . Then we are lead to consider strings in the background of D2 branes. The large 



iV limit of this theory was considered in |12[] , where it was observed that the supergravity 



description is valid only in some region of the solution. Therefore the analysis of the Wilson 
loops will also be a bit more involved. We will find that we can calculate the Wilson loops 
from supergravity only when the size of the loop is not too small. This is just related to 
the fact that for small distances we can trust the perturbative super- Yang-Mills theory. 
The physical result is quite different when the Wilson loop is large. If we consider a string 
worldsheet, embeded in the p-brane solutions studied in [1^] in a configuration appropriate 



for studing a the quark-antiquark forces we find that we have to minimize the action 

S = J dx^(d x U) 2 + U 5 /R 5 (5.5) 



where R 5 = &Tx 2 g Y M^ • We obtain solutions very similar to (O), which lead to the 
potential 

2*^(4)5/3 (g* M N)W _ r(f)t7 
3V3r(A)5/3 L 2/3 ^ r( ^) 

between quarks and antiquarks. Uq is the minimum value of U. Now we perform the 
analysis of when we can trust ( |5.6| ). Let us first consider the large U region. According 
to |12| we can trust supergravity for U <C <7y M iV. The solutions to (|5.5| ) consist of 
string worldsheets going all the way to U = oo. However the large U behaviour of the 
solution matches that of the infinitely massive W-boson. So we will require the solution 
at U ~ 9ym-N to be very similar to that of the W-boson, i.e. we require x — L/2 <C L. 
This implies that L ^> l/(gy M N). If the distance between the quarks was much smaller 
than the above bound then we can apply perturbative Yang-Mills and we would obtain a 
potential proportional to V ~ #y M iVlog(L<7y M iV). We see that these answers match up 
to a numerical coefficient with ( |5.6| ) when both calculations break down at L ~ l/(<7y M iV). 

Now we need to see if we can trust the behaviour of the solution at small U, which 
corresponds to large distances. At small U we expect that the worldsheet of the string 
turns into an M-twobrane wrapped along the eleventh direction. If Uq 3> 9ym then we can 
trust the above results ( |5.6| ) . If Uq is smaller then we have to consider a more complicated 
situation where we have to solve the equation of the M-twobrane in the background cor- 



responding to a periodic array of M-twobranes as described in [[E_2| , this presumably could 
be done but we will not attempt to do it here. 

In summary, for the non-conformal theories one can also use the description of classical 
string worldsheets embeded in some background supergravity solution to calculate Wilson 
loops when the the size of the loop is large enough (otherwise we could use the perturbative 



description). One could consider other p-brane field theories as in [Tj] and one would 
find various constraints on the size of the Wilson loop for when one can trust the string 
worldsheet description. Of course the total size is not the only issue, we also need that the 
contour does not wiggle too much. 
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